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The chiral phase structure of the Nambu-Jona-Lasinio/Gross-Neveu model at finite 
temperature T and finite chemical potential \x is investigated using (Wilsonian) Non- 
Perturbative Renormalization Group (NPRG). In the large N c limit, the solutions of 
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NPRG with various cutoff schemes are shown. For a sufficiently large ultra-violet 
cutoff, NPRG results coincide with those of Schwinger-Dyson equation and have little 
cutoff scheme dependence. Next, to improve the approximation, we incorporate the 
mesonic fluctuations. We introduce the auxiliary fields for mesons, and then derive 
NPRG equation for finite N c . The chiral phase structure on (T, /i) plane beyond the 
leading of 1/N C expansion is investigated in the sharp cutoff limit. iV c dependence of 
chiral phase diagram is obtained. 
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1. Introduction 



Explorations of the phase structure of hot/dense QCD or its toy models are of theoretical 
interest, and it is relevant to the heavy ion experiments planned at RHIC and LHC. At finite 
temperature and density, the vacuum is expected to move to the quark gluon plasma (QGP) 
phase via the phase transition, where chiral symmetry is recovered and quarks and gluons 
are deconfined. Since the phase transition to QGP phase is a non-perturbative phenomenon, 
we need a non-perturbative analysis to understand it. Unfortunately, there are not so many 
methods to deal with such problems: lattice Monte Carlo simulation, Schwinger- Dyson equa- 
tion (SDE), 1/N expansion, e expansion. Non-perturbative renormalization group (NPRG) 
is also one of the methods for such a purpose. The effectiveness of the NPRG method in the 
non-perturbative phenomena has been investigated by many authors. 

Non-perturbative renormalization group equations describe the response from the change 
of infra-red momentum cutoff A and can be written down exactly. They are the functional 
differential equations for the Wilsonian effective action in which the quantum correction 
from the high energy modes {p > A) are already incorporated. In the practical analysis, we 
approximate the theory space, the functional space of the effective action, and project the 
renormalization group equation (RGE) onto this sub-space. By enlarging this sub-space, 
we can improve the approximations systematically. In some cases, the evaluated physical 
quantities converge fast under such a process. 3 ) This is an advantageous feature of the 
NPRG method compared to the asymptotic series, e.g. the perturbation theory, e expansion 
and 1/N expansion. 

In this article we investigate the chiral phase structure of Nambu- Jona-Lasinio (NJL)/Gross- 
Neveu (GN) model 4 )' 5 ) at finite temperature and finite chemical potential. The exploration 
of the phase diagram is of fundamental interest. If the analysis is extended to QCD, it is 
of use for the early universe, the astrophysics of neutron stars and the physics of heavy 
ion collisions. We employ NPRG method for the analyses. It is worth while to examine 
the applicability of NPRG method, because there does not exist so many other tools for 
non-perturbative analyses and they do not necessarily work well in any situation. 

In lattice Monte Carlo calculation, much exploration has been made in the system at 
finite temperature, 6 ) while the behavior at finite density is much less understood. The 
non-vanishing chemical potential /i makes fermion determinant to be a complex number, 
and therefore straightforward Monte Carlo methods can not be applied. At present, two 
known candidates avoiding this difficulty, i.e. Glasgow algorithm 7 ) and imaginary chemical 
potential /i = iv method, 8 ) require much larger computer resource but unfortunately do not 
bring any definite results. Most of the efforts have been done using SDE, but it is difficult 
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to improve the approximation systematically. 

This paper is organized as follows. In § 2 we derive the evolution equation which is 
one of the non-perturbative renormalization group equations and explain local potential 
approximation (LPA). 9 ) In §3 we discuss the chiral phase structure in the large N c limit, 
and compare with the results from SDE for the fermion mass function S(q). Due to the 
formal equivalence of two methods in the large N c limit, our large N c result should coincide 
with that of SDE. We show consistent results can be obtained in the framework of NPRG 
method. In § 4, we investigate the phase structure beyond the large N c approximation. Non- 
perturbative renormalization group method can approximately incorporate the higher order 
diagrams in 1/N C expansion within local potential approximation. The N c dependence of the 
chiral phase diagram will be presented there. § 5 is devoted to the summary and discussions. 

§2. Evolution equation and local potential approximation 

There are three formulations of the non-perturbative renormalization group, the Wegner- 
Houghton equation, the Polichinski equation and the evolution equation. 10 ^ " They are 
the continuous version of the block spin transformation written in the momentum space and 
describe the response of lowering the infra-red momentum cutoff A. We can find the exact 
form of renormalization group equation (RGE) for the Wilsonian effective action and/or the 
effective average action. The latter one is the one particle irreducible part of the Wilsonian 
effective action. By lowering the cutoff we have the effective action at large distance and 
incorporate the radiative corrections from the high energy modes. In this article, we employ 
the evolution equation in Ref.12) and apply to the Nambu-Jona-Lasinio (NJL) /Gross- Neveu 
(GN) model 4 )' 5 ) at finite temperature T and chemical potential /j,. 

The generating functional of connected Green functions is 



W A [ri, V}=\nj D$Dip exp {-,S c f ut A $, ^ ~ <W?/S 1>] + fj • *l> ~ $ • v} , (2-1) 



where S, 



cut 



A is given by, 



Scut A ["0,-0] 



L 



1/fcT 




id,A)ip. *) 



(2-2) 



Here A f is a cutoff operator and has the property, 




for p > A 
for p < A. 



(2-3) 



*) In the following, we choose units such that Boltzmann constant k = 1. 
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At finite temperature, p is quantized to the Matsubara frequency, 



u)f iTl = (2n + 1)ttT or uj h ^ n = 2miT, (2-4) 

which is for fermions and bosons, respectively. Note that, since Sl ut A preserves the chiral 
symmetry, the effective action also respects it. We choose the cutoff function C as, 

C-\p/A)= f2{ fJ ( A \ v with f(p/A)=ex P {-a(p/Af% (2-5) 
l-f 2 (p/A) 

Since the parameter a can be absorbed in a redefinition of the cutoff A, we fix it such 
that a = 0.3 in this paper. The parameter b is the cutoff scheme parameter. *' Taking 
the derivative of Eq.(2-1) with respect to A and performing the Legendre transformation 
ip] + ip ■ Aj ip — fj ■ ip — ip • rj — Wa]?], T]), we get the evolution equation for the effective 
average action, 



-i f a-1 , n( 2 ) 



(2-6) 



dA' l ^ 1 

where str is super-trace which involves momentum (or coordinate) integration, Matsubara 
summation, spinor summation and color summation. r\ is a second (functional) derivative 
with respect to the fields <P = (ip T , ip) i.e. 

(^ 2) ) =^r A [*]JL. (2.7) 

V Jxy S$ x 5<P T y 

This RGE possesses the exact information about the response of the effective average action 
to the coarse graining. However, it is a functional differential equation and we can not 
solve it without approximation. As a first step of our approximation, we neglect higher 
derivative terms and keep Z factor to be unity. It is called the local potential approximation 
(LPA). 9 ) This is the leading order of the derivative expansion. 13 ) Since the NPRG preserves 

1/Nc 1/Nc 1/Nc 1/Nc 1/No 1/Nc 

+ ^rr + ~ 



(a) order 1/Nc (b) order 1/N? (c) order 1/Nc 2 

Fig. 1. The Feynman diagrams of four-fermi (3 functions incorporated in the LPA. The leading 
contribution in 1/iVc expansion corresponds to the first diagram (a). Lines denote the flow of 
color indices. 

the homogeneous global symmetries, the effective average action r A also respect it. So the 
operator space is restricted to that of the chiral invariants. For example, the independent 
*) In § 3, we employ this smooth cutoff regularization to see cutoff scheme b (in)dependence. 





4 



chiral invariant four-fermi operators are given as; 0\ = {ipip) 2 + (ipi'-f^ip) 2 , O2 = (ip'-fiip) 2 + 
(?/>7j75-?/>) 2 and O3 = (^70 ip) 2 + (^7075^) 2 Here color indices are omitted and i is the 
spatial index running over i = 1, 2, 3. In the LPA, the effective average action is 

r A $M = fJ T dr J d d ~ l x - 1^ - ]T jfi-Oi + ■ ■ j , (2-8) 

where N c is the number of colors. In the LPA, the four-fermi operators do not receive 
any corrections from the multi-fermi operators other than the four-fermi operators. The 
Feynman diagrams corresponding to the (3 functions of the four-fermi coupling constants are 
drown in Fig. 1. In the large N c limit, the Feynman diagrams (b) and (c) in Fig. 1 don't 
contribute to our (5 functions. 



§3. Phase structure in the large N c limit 

In this section, we explore the phase structure of the Gross-Neveu model 5 ) at finite 
temperature and chemical potential in the large N c limit. We attempt to apply the NPRG 
method to GN model at T 7^ and fi ^ 0, and show NPRG reproduces consistent results 
with SDE. 14 ) In the large N c limit, only the first diagram (a) in Fig. 1 contributes to the 
(5 function of the four-fermi couplings. The j3 function of the scalar four-fermi coupling 
C7 S = G\ is a function of C7 S alone, 

^C s = -(d - 2)C7 S + 2C7 s 2 /(a, b ;f,fi), (3-1) 

where t is the cutoff scale parameter i.e. A = A exp(—t). A profile of the threshold function 
I (a, b; T, p) is given in Appendix A. The characters with hat are the dimensionless coupling 
constant, e.g. G s = G s /A d ' 2 . The first term of the right hand side (RHS) of Eq.(3-1) 
corresponds to the canonical scaling and the second one to the radiative correction. At zero 
temperature and zero chemical potential, we find the two-phase structure by solving the 
RGE numerically. The RG flow diagram is shown in Fig. 2. There are two phases divided 
by critical coupling, the strong coupling phase and the weak coupling phase. In the strong 
coupling phase, the four-fermi coupling constant blows up to the infinity at finite scale t. By 
evaluating the effective potential of meson fields, we can recognize that the chiral symmetry 
is spontaneously broken in this phase. 2 - ) In the weak coupling phase, the four-fermi coupling 
goes to zero and chiral symmetry is not broken. 

Let us discuss hot and dense matter. The broken chiral symmetry at zero temperature and 
zero chemical potential is restored at some critical temperature T c and/or critical chemical 
*) We consider only the case with a single flavor. 
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Fig. 2. RG flow diagram for four- 
fermi coupling G$ of three dimen- 
sional GN model in the large N c 
limit at zero temperature and zero 
chemical potential. The cutoff 
scheme is b = 2.4. 
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Fig. 3. RG flow diagram for four- 
fermi coupling G$ of three dimen- 
sional GN model in the large iV c 
limit at finite temperature and 
zero chemical potential. The cut- 
off scheme is b = 2.4. 



potential /i c . We show the RG flow of the four-fermi coupling constant at finite temperature 
in Fig. 3. There is the critical temperature, below which the four-fermi coupling constant 
blows up to the infinity. On the other hand, above the critical temperature the four-fermi 
coupling turns to decrease exponentially. 

The temperature/chemical potential dependence of various quantities can be found by 
solving RGE with the same initial condition as that of T — /j, — 0. There are the critical 
temperature/chemical potential, above which the four-fermi coupling tends to zero, so that 
chiral symmetry is restored. We can estimate the critical temperature and the critical 
chemical potential by solving the RG flow equation for the scalar four-fermi coupling constant 
with some fixed initial condition, or equivalently bare coupling constant. The initial condition 
should be given at A — > oo. However, then we must calculate contributions of infinite 
number of Matsubara modes since T/A = T goes to zero. *-* In the practical analysis, 
we solve the RG flow equations with the common initial condition at the sufficiently large 
but finite ultra-violet cutoff Aq. If Aq is sufficiently large compared with T and fi, the 
solutions will reach the scaling region where the renormalized information i.e. Aq — > oo 
limit is obtained. This corresponds to tuning the bare four-fermi coupling constant to the 
critical one : G s \a=a = G* + 5G S , 5G S — > 0, where G* is the critical coupling constant at 

*) As seen from the explicit expression of threshold function in Appendix A, the contributions of high 
Matsubara modes are suppressed exponentially for a finite T. So at some high Matsubara modes, they 
become to make no contribution within the accuracy of numerical computation. In T — > limit, however, 
infinite number of Matsubara modes contribute. 
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T — /j, — 0. For sufficiently large Aq (or sufficiently small 5G S ), the critical temperature 
and critical chemical potential proportionally depend on the common factor 5G S d \ 
i.e. Tc(SG s ) oc SG S d \ /j,c(SG s ) oc 5G S d \ *-* Thus the ratio of these quantities 
is independent of the initial value of the four-fermi coupling for the scaling region. The 
practical problem in the calculations is how large Aq is needed to reach the scaling region. 
In Figs. 4 and 5, temperature and chemical potential dependence of threshold function / is 
shown. At small dimensionless temperature T/A < 0.03, threshold functions seem almost 
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Cutoff Scheme b=5.4 




T/A 



Fig. 4. Temperature dependence of 
the threshold function I at \i = 
with various cutoff schemes in 
three dimensional GN model. 




T/A 



Fig. 5. Temperature dependence of 
the threshold function 1 at fi 7^ 
in three dimensional GN model. 
The cutoff scheme is b = 2.4. 



constant **) and have almost no temperature dependence. So we expect to reach the scaling 

region if we set the initial condition at A > 33 T. However, in practice, we must check 

whether the physical quantities scale correctly or not, since it may depend on the quantities 

to calculate. We perform the scaling check by calculating a ratio of critical temperature (or 

chemical potential) and m . Here m is dynamical fermion mass at T — /i — 0. 

We show the results in Fig. 6. If we employ the ultra-violet cutoff A = 333 T c „, the 

critical points (T c /m , fic/m ) almost coincide with those of SDE. Here T Co is a critical 

temperature at /j, — 0. There is little cutoff scheme dependence using this ultra-violet cutoff 

A . With the ultra-violet cutoff A = 333 T c these quantities reach the scaling region. 

On the other hand, if we choose the ultra-violet cutoff A = 33 T c , the critical points 

(T c /m , pLc/mo) deviates from those of SDE a little. A = 33 T c is not sufficient to obtain a 

scaled quantities. We can calculate scaled critical temperature/chemical potential by setting 

the initial condition at A > 333 T c „■ We also plot the results in the sharp cutoff case i.e. 

*) Other physical quantities, such as the fermion effective mass mo at T = fj, = also proportionally 
depend on the factor SGg 1 ^ 2 ^^. For the detailed explanation, see Appendix B. 
**) They are the values at T = \i = 0. 
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b = oo in Fig. 6. *•* 

It should be noted that above analyses are based on the flow equation for 1/G S rather 
than Eq.(3T) for a technical reason. In the large N c limit, 1/G S coincides with the mass 
squared of the meson fields. Multiplying — 1/Gf to RGE (3-1), it can be rewritten as 



(3-2) 
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Two phases, chiral broken phase and symmetric phase on (T, ji) plane can be distinguished 
by the behavior of the solution of eq.(3-2). If 1/G S goes to the negative value at suffi- 
ciently large t, it is the strong coupling (broken) phase. If l/G s tends to the positive 
infinity, it is the weak coupling (symmetric) phase. Although Eq.(3-2) and Eq.(3T) are 
essentially identical, there emerge difference 
on the analyses. If chemical potential j2 is 
large enough, the threshold function / can 
take negative value. So even if 1/G S become 
negative, l/G s may turn back to the pos- 
itive value due to the negative /. Thus it 
should be regarded as weak coupling phase. 
With RGE Eq.(3-1), this behavior can not 
be detected because once G s blow up to the 
infinity, then by mistake we may identify it 
as the strong coupling phase. Indeed, for a 
large /x, it gives an inaccurate (about 10%) 
result. The analysis using 1/G S criticality is 
a more favorable method and we may say, it 
is a sort of "environmentally friendly renor- 
malization group" . 15 -* 

In Ref. 14), the chiral phase structure of 
GN model in 2 < d < 4 dimensions were 
calculated using SDE. The effective poten- 
tials are calculated in the leading order of 

the 1/N C expansion. It is known that in 2 < d < 3 chiral phase transition is of first order at 
low temperature and large chemical potential, and of second order in other region of critical 
line. In 3 < d < 4, the phase transition is of second order all along the critical line. 

Using NPRG we can also derive consistent phase boundaries in any space-time dimen- 
*) The detailed explanation in the sharp cutoff case is given in § 4. 



Fig. 6. Chiral phase structure of three dimen- 
sional GN model at finite temperature and 
chemical potential in the large N c limit. 
Temperature and chemical potential are 
normalized by dynamical fermion mass niQ 
at T = n = 0. The dependence of both 
the cutoff scheme and ultra-violet cutoff A 
in the case of NPRG method Eq. (3.2) and 
also the results of SDE are shown. The 
phase transition is of second order all along 
the critical line. 
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sions. For example, the chiral phase structure of two dimensional GN model *) is shown in 
Fig. 7. In two dimensions, there is the first order critical line as well as the second order 
one. The phase boundary can be found only by exploring the effective potential V of the 
collective coordinate a, which is introduced in the next section. The second order critical 
point is where d 2 V/da 2 \ a= o = is satisfied, while the first order critical point is where the 
value of the effective potential at the second minimum is equal to that at the first minimum 
at the origin. The stars in Fig. 7 stands for the first order chiral phase boundary and it 
terminates at the tricritical point. The filled circles and the open circles are where the signs 
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Fig. 7. Chiral phase structure of two 
dimensional GN model at finite 
temperature and chemical poten- 
tial in the large Nq limit. Temper- 
ature and chemical potential are 
normalized by dynamical fermion 
mass mo at T = /i = 0. The sym- 
bols are the results using NPRG 
method. Here cutoff scheme is 
b = 1 .4. Ultra-violet cutoff Aq is 
333 T c I,. The lines are the results 
using SDE. Phase boundary con- 
sists of the first order critical line 
and the second order one. The 
dashed line and the open circles 
have nothing to do with the phase 
transition. 
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Fig. 8. RG history of the effec- 
tive potential of two dimensional 
GN model on the first order 
phase boundary in the large iVc 
limit. We choose T/mo = 
0.113, n/mo = 0.682. The effec- 
tive potential and the field ex- 
pectation value are normalized by 
the dynamical fermion mass mo at 
T = /j, = 0. Cutoff scheme here is 
b = 1.4. Ultra-violet cutoff Aq is 
333 Tq o • 



1.4 



of meson masses at the origin < a >= are changed and they correspond to the 1/G S 



*) In two dimensions continuous symmetry can not be spontaneously broken. 16 ) So we consider the 
phase transition with a discrete chiral symmetry if) — > 75^ in two dimensions. 
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critical line. The filled circles are identical to the second order chiral phase boundary, but 
the open circles have nothing to do with the phase transition. The effective potential of two 
dimensional GN model has a gap structure on the first order phase boundary. An example 
of the RG evolution of the effective potential at a first order critical point is shown in Fig. 8. 



§4. Phase structure beyond the 1/N C leading 

In the previous section, we explored chiral phase structure in (T, /x) plane in the large 
N c limit. In this section, we improve the approximation beyond the 1/N C leading in four 
dimensional NJL model. Since we are interested not only in the phase boundary but also in 
the order of phase transition, the CJT effective potential should be investigated. Let us 
start from the definition of the partition function, 

Z A [J] = J D4>D4> exp {-S c f ut A - S bare + fj ■ i/> - $ ■ rj + ■ $i/> + K ■ ^7s^} , (4-1) 

where S { cut A , S^are and J are the cutoff action defined in Eq.(2-2), the bare action and the 
sources J = {i], fj, E a , S w }, respectively. Here the bare action is taken to be, 

Share = / d 4 X - " ^ [(W + (fal^f] } ■ (4-2) 

We introduce the auxiliary fields a, ir corresponding to the composite operators ipip, ipij^ip 
by the following Gaussian integral, 

Z A = f DaDit e~ Stt A ~^ x {h^-^-P^TH^-^-P^T} ^ (4 . 3) 
where S^ut A is the cutoff action of meson fields and given by, 

^cut A = Jd A x l - (d,a C-\-id/A) d,a + d,n C~\-id/A) d,n) . (44) 

Here C™ 1 is the cutoff function defined in Eq.(2-5). Note that, since Z A depends on the fields 
■0, ip as well as on the sources £ a , £ n , the partition function is deformed by the insertion of 
Z A except at A — 0. Inserting Z A into the path-integral of Eq.(4-1), the partition function 
becomes 

Z' A [J\ = j DipDipDaDn e~ s ^ *" s «t ^s^+^-^+p-^e^+p-^^ 

xe " / dAx {I^M^ 2 +^ 2 )-^(^+^75)^+ip[W) 2 +(^75V') 2 ] + 2^(^ 2 +^ 2 )} 
-fD*D<, e-«- (4-5) 
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where is given by, 

Sbarc = ^barc ~\ j ^{ fa + P~* Z„) [(P" 1 + C^U)- 1 - P] fa + P^E, 

+ (i>i lb i> + P' l E^) [(P- 1 + C^uy 1 - P] (^ 75 V + P" 1 ^) |. (4-6) 

In the limit A ^ 0, Z' A = Z A except for the field independent constant. We modify the gen- 
erating functional W A [J] = \nZ' A [J] as W' A = W A + 1/ '(2P)-(Z% + Z%). Physically meaningful 
quantities are not affected by this modification. 18 ) Indeed, no vacuum expectation values is 
changed by adding an arbitrary polynomial of the sources to the generating functional, e.g. 
in our case, 



dW A 



dW' A 



dS„ 



- P- 1 ^ 



dW' A 



So 



d£„ 



(4-7) 



We rescale the sources as J' = {77, 77, E' a , E' n } = {77, r/, P 1 £ cn P 1 E n }. The initial condition 
of RGE at the ultra-violet cutoff A — A can be given after the Legendre transformation, 
r A C3T [ip, if>, a, tt] = f] ■ ij> - $ ■ r] + E' a ■ a + 27; • 71 - W' A [J'\. We have 



pci-i _ _|_ ^ _|_ ^ _|_ 



±j d^p- 1 {(a - p^y + (tt-p^^) 2 } 



= S£t\ + / - V*70 - (7 - « 75 )^ + ^(a 2 + TT 2 )} , (4-8) 

where S^ A = S { cut A + 5 c b ut A . On the last line of Eq.(4-8), we took P = G s /N c . By 
subtracting the cutoff action we define the effective averaging action as : r A [ip, a, tt](— 
J d 4 xC A ) = r A iT — A . In this paper, we approximate it as follows : 



C A = i>(Zf id - ifx)j ip + Zf^idai^j + -^N c Z m ((d cr) 2 + (<9 7r) 2 ) 
+ lN c Z b ((d t a) 2 + (<9,7r) 2 ) + y$(<r + i 75 7r)^ + N c V(a 2 + tt 2 ). 



(4-9) 



There are a yukawa interaction and an effective potential V which is a function of a 2 + n 2 . 
Here, by a lack of Lorentz symmetry, the wave-function renormalization factors of temporal 
derivative and those of spatial one are different each other. 

The RG flow equation for the effective average action P A can be derived as done in 
§2. In our approximation, the RGE becomes a partial differential equation for the effective 
potential V(p) where p = 1/2 • (a 2 + rr 2 ). If we employ the smooth cutoff function C(q/A) 
as done in the previous section, we have to perform the momentum integration with respect 
to the spatial momenta numerically for each Matsubara modes. It is very time-consuming 
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and not practical. As investigated in the previous section, the results do not depend on 
the cutoff scheme for a sufficiently large ultra-violet cutoff ylo at least in the large N c limit. 
Therefore it is more promising way to take the sharp cutoff limit of the RG flow equation. 19 -* 
After doing so, we need no momentum integration and therefore the computational time is 
reduced drastically. We can not do Taylor expansion in momentum components p^ because 
the sharp cutoff induces non-analyticity at the origin of the momentum space. 19 **' 10 - ) Instead 
we expand one particle irreducible generating functional F A in momentum scale \p\ = ^Jp^p^- 
*) 

In general we may express the cutoff function C(q/A) as C~ 1 (q/A) = l/9 £ (\q\, A) — 1 in 
terms of # £ (|g|,/l) which is a smooth regularization of Heaviside 6 function of width e. In 
the sharp cutoff limit (e — * 0), e (|<7|, A) reduces to 9(\q\ — A) and the evolution equation for 
the effective average action becomes 

(4-10) 



-^r A [ijj^,a,ir] = ~str 



6(\q\ - A) ~ (2) , r A(2)\-i 

———r A [i + cr A ) 



where we separate the field independent full inverse propagator j{p, A) from the two point 
function 



(ri 2) [^^,a,n]) ppi = -y(p,A) (2vr) 4 5(p + p') + (jf ^ a, tt])^ , (4-11) 

so that /1 2) [0, 0, 0, 0] = 0. 19) In Eq. (4-10), G(p, A) is the infra-red cutoff propagator, 

G(p, A) = hm - u -i -— = e ^ P )~f ) . (4-12) 

e^oC-\p/A)+ 1 (p,A) -y(p,A) 

As mentioned above, in the sharp cutoff limit the origin of momentum space is not analytic 
and Taylor expansion in p^ breaks down. Alternatively, one should expand in momentum 
scale. The infra-red cutoff function is expanded in terms of absolute value of external mo- 
menta \p\ as 

oqp + q \- A) = e(q.p+ \ p \/2) = e( q ■ p) + f; -1 a*™- 1 ^ ■ p). (443) 

m =i m. \ Z J 

Here p^ is a unit vector parallel to p^. Integrating RHS of Eq. (4-10) with respect to 
the internal momenta q, one can expand it in terms of the momentum scale \p\. In finite 
temperature case, the momentum integral in Eq. (4-10) turns to the angular average over 
three dimensional momenta which are restricted to q 2 = A 2 — uo 2 b or ^ n for each Matsubara 
mode. Using the formula of angular integral, 

/ d 3 q 8 (|q| - ^ - cof b or /)jn ) /(q • p) = 2n(A 2 - uf h or /}j J £ dz f(z), (4-14) 



*) The Mth order approximation is to drop all terms beyond 0(\p\ M ) and the lowest order approximation 
coincides with the local potential approximation of the Wegner-Houghton equation. 19 ) 
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we can integrate with respect to the internal momenta q and expand RHS of RGE 
(4-10) in terms of the external momenta |p|. The flow equation for the effective potential 
V(p) can be derived. We have 



d_ 

dt 



V = (3 f + (3 b , 



(4-15) 



where flf is the contribution from fermion loop graphs and [3 b is that from boson loop graphs. 
They are expressed as; 



T 



7T 

T 



1 + 



'y p 



An 2 N c n 
T 



z2'(b,n In 



ZjA 2 + (Zj - Zj) u} tn -p 2 + 2i Zfo u f , n p 
Z b A 2 + (Z b0 - Z b ) u 2 + V' 



4ir 2 N ( 



]T'C M In \z b A 2 + {Z m - Z b ) u{ n + V + 2pV" 



C n 



T 



An 2 N, 



■J2'Cb,n In \z b A 2 + {Z m - Z b ) uol n + V 



C n 



p=0 



(4-16) 



where (( b or f^ n = ^ A 2 — u 2 b or ^ ~ • A 2 and the prime ' operating on the effective potential 
V denotes the derivative with respect to p. J2' n is the summation about n with the condition 
fb or f),n — A 2 - For the wave function renormalization factor Z b , we have *^ 



dt b vr^ {Z)A 2 + (Z% - Zj) W J >B - p 2 - 2iZ fQ p u f , n ) 2 
(I 1 Z}(A*-u} >n ) \ 



Cf,nZj y 2 



(4-17) 



\,2 3 Z 2 A 2 + (Z% - Zj) ul n -p 2 - 2tZ f0 p u ftn J ' 

In this paper, we also approximate Z b Q as Z b o = Z b . In the large N c limit, the boson 

fluctuations disappeared, i.e., f3 b — due to a factor 1/N C . However for finite iV c , we must 

incorporate the boson loop contributions, f3 b . This is a significant difference from the large 

7V C limit. In general the fermion contribution makes the effective potential evolve downward 

as cutoff is lowered, while the boson contribution lifts the effective potential upward. **) The 

latter is main additional effect for finite N c on the RG evolution of the effective potential 

and hence also on the chiral phase structure on (T, p) plane. We shall ignore the corrections 

to fermion's Z factors Z/, Zj and to yukawa coupling y***^ as the first step toward the 

*) In momentum scale expansion of sharp cutoff effective average action, there emerges a kinetic term 
proportional to |p|. The non-analytic dependence on momenta at p = lead to a non-locality in position 
space. Since the resulting physics in the limit A — > should not suffer from this non-locality, it should be 
absorbed by a certain counter term. Here wc simply neglect it in our approximation. 
**) This is understood from the signs of /?/ and (3 b - 

***) Here, we do not renormalize the boson fields a, 7r to make their Z factor unity, i.e. Z b = 1, since at 
ultra-violet cutoff, it should vanish. Therefore in an ordinary sense, our physical yukawa coupling runs due 
to bosonic wave-function renormalization. 
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exploration of the phase structure at finite temperature and finite chemical potential beyond 
the 1/N C leading. *) 

RGE for the effective potential is non- 
linear partial differential equation. Here we 
do not attempt to expand it in powers of the 
fields as done in Ref.3), since if the cubic 
term of the field i.e. p 3//2 appears in the effec- 
tive potential V(p), 21 ^ then the naive poly- 
nomial approximation of it can not work. 
We discretize the RG flow equation in t and 
p directions and solve it numerically. We 
apply the extended Crank-Nicholson method 
22), 23) w ] 1 i c ] 1 j s effective for a non-linear par- 
tial differential equation. 

By evaluating the minimum of the effec- 
tive potential, we can conclude whether it is 
a broken phase or a symmetric one. If there 
is a non-trivial minimum at nonzero value of 
a, it is a broken phase. If the temperature (or the chemical potential) is greater than the 
critical one, the effective potential has a minimum at the origin, i.e. the chiral symmetry 
is recovered. If the temperature (or the chemical potential) is below the critical one, the 
effective potential has a nontrivial minimum. More precisely, in the broken phase, the effec- 
tive potential evolves toward the so-called 'convex' one. 25 - ) ' 26 - ) Some results of the effective 
potentials are shown in Fig. 9. 

Table I. The critical temperature T c q and the chemical potential [i c q of the four dimensional NJL 
model. 





10 


20 


40 


80 


oo 


Tco/mo 


1.07 


1.21 


1.30 


1.35 


1.38 




2.58 


2.70 


2.81 


2.83 


2.83 


A*c o/Tc o 


2.41 


2.23 


2.16 


2.10 


2.05 



Taking into account of the discussions in the previous section, we take an initial condition 
of RG flow equations at A = 333 T Co . Indeed, with this condition, the critical line almost 
coincided with that from SDE in three dimensional GN model for N c = oo (See Fig. 6). 
*) A similar approximation is performed in Ref. 20). 




3.0 



o/m. 



Fig. 9. The effective potential of four dimen- 
sional NJL model at finite temperature. 
Here, iV"c = 10. The effective potential 
and the field expectation value are normal- 
ized by the dynamical fermion mass mo at 
T = n = 0. 
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0.0 0.5 1.0 1.5 

T/m 

Fig. 10. N c dependence of the chiral phase structure of the four dimensional NJL model at finite 
temperature and chemical potential. In each N c , the phase transition is of second order or 
very weak first order all along the critical line. The temperature and the chemical potential 
are normalized by the dynamical fermion mass mo at T = /j, = 0. 

The chiral phase diagram on (T, /i) plane is shown in Fig. 10. Here, the temperature and 
the chemical potential are normalized by the dynamical fermion mass m = y < cr > at zero 
temperature and zero chemical potential. In each N c case, the phase transition is found to 
be of second order or very weak first order all along the critical line from th shape of the 
effective potentials at the critical points. For example, an effective potential at a critical 
point for N c = 10 is shown in Fig. 9. N c dependence of the critical line on (T, /i) plane is not 
small. The critical temperature and the critical chemical potential become small in unit of 
m if N c is lowered. This is due to the boson fluctuations existing for finite N c . Intuitively, 
we can understand it as follows. The fermionic negative corrections to the effective 
potential are suppressed by a high temperature and/or a large chemical potential, but the 
bosonic positive ones ((3b) are less suppressed due to the existence of Matsubara zero mode. 
Hence if the bosonic part becomes large, that is, the case of small N c , the chiral symmetry 
is restored at lower temperature (smaller chemical potential) in comparison with the large 
N c case. N c dependence of the critical temperature at /i — (T c „) as well as that of the 
critical chemical potential at T = (// c ) are shown in Table I. *) 



*) The numerical calculations for N c < 10 need much finer mesh in t direction because /?(, becomes 
large. This will be a subject in the forthcoming paper. 
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§5. Summary and discussion 

We performed a non-perturbative analysis by the (Wilsonian) non-perturbative renormal- 
ization group in NJL/GN model at finite temperature and chemical potential and explored 
chiral phase structure on (T, /i) plane. 

First we explored the phase structure of the GN model at finite temperature and chemical 
potential in the large N c limit. We showed that consistent results with those of SDE could be 
obtained in the framework of NPRG method. We investigated how large ultra-violet cutoff A 
is needed to obtain scaled quantities. We could calculate scaled critical temperature/chemical 
potential by setting the initial condition at A = 333 T c and they had little cutoff scheme 
dependence. 

Second we improved the approximation beyond the 1/N C leading in four dimensional NJL 
model. For finite N c , the bosonic fluctuations have to be incorporated. In §4, we introduced 
the mesonic auxiliary fields a, ir, and derived the flow equation for the effective potential of 
these fields V(a 2 + ir 2 ). RG flow equation is then a non-linear partial differential equation 
for this potential. Taking account of the possibility of the first order phase transition like 
Fig. 8 and the cubic term p 3 / 2 , 21 ) we solved the partial differential equation for V(a 2 + it 2 ) 
without a polynomial expansion employed in Ref.3). For several N c , we obtained the chiral 
phase structure on (T, ji) plane. For N c = 10 ~ oo, the phase transition is second order or 
very weak first order all along the critical line. The critical temperatures/chemical potentials 
depend on N c largely, and become small as N c is decreased. For N c = 10, T c is 78% of 
that for N c = oo. 

Generalization of this analysis to more realistic models, QCD seems to be straightforward, 
except for the treatment of the gauge invariance. By introducing the instanton induced 
multi-fermi operator, 24 ) NPRG method can be also applied to the color-superconductor. 27 ) 
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Appendix A 

Threshold functions 

In this appendix we show the explicit form of the threshold functions I(a, b ;T,jT) ap- 
peared in Eq.(3-1). In d dimensions, the threshold function at zero temperature and zero 
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chemical potential is 

I r°° 

/(a ' 6 ;0 ' 0) = ¥=^rm L dq ^ 2{d/2+2)a b - / 2 (?))/ 2 (?)- (A-i) 

For finite temperature and/or chemical potential, the threshold function can be obtained by 
the replacement / dq d /(2n) d ^ T dq^ 1 / (2n) d - 1 and q -^Q f , n -ip = (2n + l)nf - ip, 

•(1 - / 2 (g))/ 2 (g)T ^ n T + g2 v ( A '2) 

fe + (? + f}H n + («-/i) 2 } 

f -> and /2 -> limit of Eq.(A-2) coincides with Eq.(A-l). 

Appendix B 

Renormalization and the continuum limit atT,/i^0 



Let us discuss renormalization and the continuum limit at finite temperature and finite 
chemical potential in this appendix. In the large N c limit, the RG (3 functions are given by, 

^G s = -(d-2)G s + I(f,fi)Gi (B-l) 
|f = f, (B-2) 

Jt n = ^ (B ' 3) 

Fixed points of the above RG equations at T — p — are the Gaussian fixed point G s = 
and Gross-Neveu point G s = (d - 2)/I(f = 0,p = 0) = G*. Note that, f = oo and/or 
p — oo are also fixed points. If T = oo and/or p — oo then only Gaussian fixed point 
can be found for the finite G s , since I(T,p) vanishes. The difference from an ordinary zero 
temperature and zero chemical potential field theory is T, p dependence of the threshold 
function /. As will be explained below, it does not affect the 'renormalization'. 

We linearize the RG equations around the GN point. *) Letting G s = G* + 5G S , we have, 

j5G s = -{d - 2) 5G S + 2/ Gl 5G s + (f J T + pig) Gf (B-4) 
|t = T, (B-5) 
— p = p, (B-6) 



*) For the infinitesimal neighborhood of the fixed point, this linearization is valid. 
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where I = I(f = 0, fi = 0), J T = (dl/df)(f = 0,/t = 0) and iff = (dl/dp,)(f = 0,/t = 0). 
By Lorentz symmetry at T = /t = 0, Iq = 0. We also find /J = 0, since the difference 
between J(T, 0) and Jq is 0(T 2 ). *) Consequently, we find, 



/ 


5G S 


\ 


( ~{d 


- 2) + 2I G$ 


\ 


' 5G S 


\ 




f 








1 


f 




V 






V 





1 / 


v £ 





The eigenvalues are found to be v\ — d — 2, i/ 2 — v% — 1. Here we used G* = (<i — 2)/ Iq. We 
derive the relations between the bare coupling and the physical quantities using an ordinary 
procedure. Let us focus on the infinitesimal neighborhood of GN point, i.e. we consider 
the continuum limit of the theory Aq — > oo. The critical temperature T c is completely 
determined by mutually independent three variables: the bare four-fermi coupling 5G S , bare 
chemical potential fio and the ultra-violet cutoff Aq. The extra jlo dependence of the critical 
temperature i.e. T C (5G S , Aq; /to) does not affect to the following considerations. Now it is 
convenient to think that the critical temperature depends on the 'bare' chemical potential /to 
through the fixed as well as RG invariant dimensionless ratio JIq/Tq. As everyone knows, the 
RG preserves the physical quantities. Therefore, we can realize by solving the RG equations 
from Aq to A — \Aq that the critical temperature satisfies, 

T C (5G S , A ; fo/T ) = T C (\-^5G S , XA ; /t /f ), (B-8) 

where we used the solution, 

SG S = Ce {d - 2)t = C(A /A) d - 2 . (B-9) 

Here C is the 'bare' 5G S (A ). The relation (B-8) can be rewritten by the dimensionless 
critical temperature f c (5G s ; p, /f ) = T C (5G S , A ; fi /f )/A **) as 

f c (5G s ; fa/To) = AT C (A-^G S ; /to/To). (B-10) 

Hence, 

f e (6G a ;fio/f )<x6& B 'W. (B-ll) 
If one starts from the temperature dependent critical chemical potential, then one finds, 

fi c (SGs; fio/To) oc 5Gl /id ~ 2) . (B-12) 

*^ I(T, 0) is like an approximation of the integral I by a histogram, band width of which is 2ttT. As 
a trapezoidal rule estimate differs from an integral by an order of the square of the width, 22 ) I(T, 0) also 
differs from I by 0(f 2 ). 

**1 By dimensional analysis, T c is independent of Aq. 
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One can also find a similar relation of the fermion dynamical mass at T — [i — 0, 
mo oc 8Gl /{d - 2) . We can 'renormalize' these quantities by taking 5G s (Aq) = (M/Aop d \ 
where M is some finite reference mass scale. Letting Aq — > oo, we find the continuum limit 
of GN model at finite temperature and/or finite chemical potential. 

The above observations can be straightforwardly generalized to other models. 
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